Spin-polarized tunneling current through a ferromagnetic insulator between two 

metallic or superconducting leads 
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Using the Keldysh formalism the tunneling current through a hybrid structure where a confined 
magnetic insulator (I) is sandwiched between two non-magnetic leads is calculated. The leads can 
be either normal metals (M) or superconductors (S). Each region is modelled as a single band 
in tight-binding approximation in order to understand the formation of the tunneling current as 
clearly as possible. The tunneling process itself is simulated by a hybridization between the lead 
and insulator conduction bands. The insulator is assumed to have localized moments which can 
interact with the tunneling electrons. This is described by the Kondo Lattice Model (KLM) and 
treated within an interpolating self-energy approach. For the superconductor the mean-field BCS 
theory is used. The spin polarization of the current shows a strong dependence both on the applied 
voltage and the properties of the materials. Even for this idealized three band model there is a 
qualitative agreement with experiment. 
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I. INTRODUCTION 



The understanding of spin-polarized tunneling cur- 
rents in mesoscopic hybrid structures is a very active field 
in current researcbi^. This interest stems from the pos- 
sibility to build electronic devices which take into account 
the spin degree of freedom of the electrons ("spintron- 
ics"). Usually magnetic tunnel junctions (MTJs) consist 
of two ferromagnetic metals separated by a thin non- 
magnetic insulator. However, in this paper we consider a 
different geometric structure. The two outer leads shall 
be either normal metals (M) or superconductors (S) while 
the central region is always assumed to be a ferromag- 
netic insulator. The band splitting in the insulator due 
to the exchange interaction will lead to different tunnel 
barrier heights for spin-up and spin-down electrons, re- 
sulting in a finite spin polarization of the tunneling cur- 
rent. This is called the spin-filter effect. 
Since a MTJ is build by coupling materials with in gen- 
eral different chemical potentials, it is inherently out 
of thermal equilibrium. Therefore the usual many- 
body description has to be modified. Its extension 
to non-equilibrium phenomena is the so-called Keldysh 
formalisms. 

Pioneering work in this field has been done by Caroli et 
al£ who first presented a model Hamiltonian for calcu- 
lating the tunneling current. The model was later refined 
by Meir and Wingreen&. Heide et alZ investigated the 
relation between the current and the interlayer exchange 
coupling. A similar model to the one presented in this 
paper was used by Zeng et al4 to calculate transport 
properties in hybrid structures. They use ferromagnetic 
and superconducting leads with a non-magnetic central 
region while in our case the situation is reversed. 
Most experiments on spin filtering were done with Eu- 
ropium chalkogenide barriers and metallic contacts. Re- 
cently there also have been several experiments where 
other materials were used. Figielski et al. 9 performed 
experiments for EuS barriers embedded into a semicon- 
ducting PbS matrix. Gajek et alJ^ used BiMn03 as tun- 



nel barrier and the half- metallic oxide La2/3Sr!/ 3 03 as 
counter-electrode. 

The paper is organized as follows: in section[TT]we develop 
the theoretical model and derive a general formula for 
the tunneling current. In sections IlIII and [TVl the model 
is evaluated for the M/I/M- and M/I/S-systems, respec- 
tivly. Special emphasis is made to understand the current 
and spin polarization features in terms of the quasipar- 
ticle density of states (QDOS). In the last section fVl we 
shortly summarize the paper and give an outlook with 
possible improvements and extensions of the model. 



II. MODEL AND THEORY 

In this section we formulate a model Hamiltonian 
which will allow us to simulate the tunneling of electrons 
through an insulating region. The Hamiltonian was first 
used by Caroli et al£. In Fig. [T]the geometric setup is 
shown schematically. At the beginning, i.e. at t = — oo, 
it is assumed that there is no contact between the three 
regions. Therefore each one is in thermal equilibrium, 
characterized by its respective chemical potential. Then 
the coupling between them is turned on, resulting in an 
overall chemical potential in the whole system. This cou- 
pling will be described by perturbation theory in non- 
equilibrium. We assume the tunneling to be ballistic, 
i.e. there is no scattering of electrons involved. Hence 
the tunneling is elastic and the energy of the electrons 
is conserved. Furthermore we assume translation invari- 
ance in each region. 

The model Hamiltonian consists of three parts: 



H — Hl + H ins + Hj 



(1) 



i?L describes the leads, Hi ns the insulating region and 
Ht the tunneling between them. For now the contacts 
are assumed to be non-interacting, thus 



(2) 



M 
M=L,R 



2 



Metal L Insulator I Metal R 




FIG. 1: Tunnel junction without applied voltage for T = 0. 
The rectangles symbolize the conduction bands of each layer. 
The chemical potential is equal for all three layers and marked 
by a dotted line. Occupied states in the metals are hatched. 

Ck Jf o- ( c k JfCr ) is the annihilation (creation) operator 
which annihilates (creates) an electron with spin a and 
wave vector ]s.m in lead M where M = L for the 
left lead and M — R for the right lead. €k M = 
Af-E, i £i M 7 M e -ik "' (R *M- R iM) is the Fourier trans- 
formed hopping integral between lattice sites R^ M and 
Rj M . In the case of superconducting leads Hl has to be 
modified. This will be discussed in section ITVl in detail. 
For the derivation of the tunneling current formula, it is 
not yet necessary to specify the interaction Hint in the 
insulator. The only restriction is that it has to commute 
with the particle number operator of the leads. Since 
Hint will in general be a functional of the insulator con- 
struction operators, i.e. H mt = H int (d^ ian : rf k f( r,m), this 
is true for all usual interactions. Hence 

Hins — e k 7 ^k/CT.n^kjCT.m + -^mt (^k/ oyre ' <^ki<T,ro) (3) 

k/cr 

The construction operators in the insulator are labeled 
by the letter d. Hi ns is written in a multi-layer form for 
reasons that will be discussed later on. So dkjo-.n (d^ I(T „) 
destroys (creates) an electron with spin a and wave vec- 
tor kj in the layer n of the insulator. 
It remains the definition of the central part of the model, 
namely the tunneling Hamiltonian Ht- Tunneling is 
modelled by a hybridization between the conduction 
bands of the leads with the (empty) conduction band of 
the insulator. Therefore we choose the following ansatz: 

R T= E E (^A M J^,n+H.C.) (4) 

M k m k / a 

M=L,R n 

H.c. is the Hermitian conjugate of the first term. 
The coupling constants e^ki are parameters which de- 
termine the strength of the hybridization between the 
bands. Their value is initially unknown and it is not pos- 
sible to determine it within the model. However, they 
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FIG. 2: Simulation of the tilting of the insulator band by d 
layers with shifted band centers. They are shifted according 
to Eq. (Hp. Occupied states are hatched and the chemical 
potential is shown as a dotted line. 

can be fixed by comparison with experiment. 
Still the Hamiltonian is not complete, since the voltage 
V does not appear anywhere. It will lead to a rigid shift 
of e • V between the band centers of the leads and the 
insulator band will be tilted accordingly. The tilting of 
the insulator conduction band cannot be modelled by a 
single insulator layer, since it can only have one center 
of gravity. Therefore the insulator is replaced by sev- 
eral insulating layers^ whose band centers are shifted 
according to Fig. [2l This is described by the additional 
Hamiltonian 

H V = V M4 Ma C tMa +J2Vrd+ a nd tI a,n (5) 

iMf il<? 
M=L,R n 

The Vm and the Vp are the potentials in the leads and 
the insulating layers, respectively. Their relation to the 
applied voltage V is as follows 

V L = 

V R = V (6) 
VP - V^j (n = l,...,d) 

This choice will lead to the behavior shown in Fig. [2] 
Since Hy only changes the one-particle energies in each 
region by a fixed amount, it can be incorporated in the 
model Hamiltonian ((T|). 

It is convenient to decompose the current flowing through 
the junction into two currents, one from the left (</£) 
and one from the right lead (J^)- Since we are only 
interested in steady state currents, the symmetry relation 
J£ — —Jr holds and the total current J a can therefore 
be symmetrized: 

J° = Jl = \{Jl + Jl) = \{Jl-J a R) (7) 
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Deriving a formula for J£ will be the next task. It will 
be proportional to the change of electron number in the 
left lead: 



(8) 



The time dependence of the current is of a strictly for- 
mal nature. The commutator on the r.h.s. has to be 
evaluated to get a final expression for J£. It is basically 
identical to the lesser Green's function G^. lkL<7 n (t,t') — 

i(c^ L(7 (t')d] lirT _ n (t)). This function can be calculated 
within the Keldysh formalism by using the analytic con- 
tinuation rules by Langreth— on the corresponding non- 
equilibrium Green's function (NEGF) G k/kl/(Tin (i, i'). 
The details of the calculation are carried out in appendix 
El One gets 



J? = 



2e f f dE 

— > Re / — e 



ktkj 



ki,kj c k L kj 



(9) 



k;(T,nm 



(E)g^(E))) 



The free metal Green's functions are given by 
g< L (E) = 2-KMf L (e kL )6(E - e k J 

9U E ) = 



E-e kL - iO^ 



(10) 
(11) 



is the Fermi function of the left 



The insulator Green's functions G k a nm (E) 
■ d+ 



where fi, (E) 
lead. 

((<W,n;< JCT ,J) G£ iaMn (E) = i(d+ i!7 ^d kia , n ) can 
be calculated by using the equation-of-motion method. 
This is done in appendix [Bj At this point of the calcula- 
tion it is necessary to specify the interaction Hi nt in the 
insulator, which is assumed to consist of localized mo- 
ments. Since there are also itinerant electrons (through 
tunneling) in the insulator, the KLM is considered to be 
an appropriate model for this region. Thus 



Hi, 



2^^^^ Za ^ii,n n il<y,n + S iin d^ i _ an di I a,n) 
n 

(12) 

where J is the exchange coupling constant and Si f tU de- 
notes the spin operator of the localized moments. The 
electrons interact only at one lattice site at a time with 
the localized moments. Therefore the summation is over 
one layer index only. The interaction causes no transi- 
tions between the layers in other words. The KLM de- 
fines a complicated many-body problem and can only be 
solved approximately for most cases of interest. Here we 
use an interpolating self-energy approach (ISA)- 3 . The 
self-energy of the KLM is in general defined as 



(([H mt ,c kia Uc+ ia )) = £ k7CT (£)G kjCT (£) 



(13) 



The ISA self-energy is an exact solution of several lim- 
iting cases of the model. It fulfills the zero-bandwidth 



limit, ferromagnetic saturation, second-order perturba- 
tion theory and the high energy expansion. The approxi- 
mation lies in the assumption that it is valid for all inter- 
mediate cases, too. It was found in other works that it 
gives reliable results for the materials usually described 
by the KLM, namely manganites, magnetic semiconduc- 
tors and local-moment metal o 141516 . The ISA self-energy 
is given by the following expression 



4 1 - UG (E- \Jm a ) 



{14) 



where 



Ni^E- e kl 

kj 



(15) 



a a = 5(5 + 1) -m a (m a + 1) (16) 
m a = z a (S z ) z (T = 5^-5 ai (17) 

The ISA is an approximation for the case of low band 
occupation. Since we use it for the description of an in- 
sulator conduction band this is no problem, of course. 
Due to the complicated structure of the self-energy, espe- 
cially the Green's function in the denominator, one has 
to be careful when applying the analytic continuation 
rules. In appendix [B] a method is shown how it can be 
incorporated in the Keldysh formalism. 



III. METAL/INSULATOR/METAL-SYSTEM 

In each of the three regions of the tunnel junction there 
is only one s-like band which is described in a tight- 
binding picture. Therefore each band is completely de- 
termined by its band center, band width and the lattice 
type. We only consider simple cubic lattices in this pa- 
per. The magnetization will be used as a model param- 
eter, i.e. it is not calculated self-consistently. It can be 
approximated by a Brillouin function^: 



(Sz 



25+1 



1 



25tanh((25 + l)^) 
35 T c 



25tanh(-|) 



5 + 1 T 



(18) 



(19) 



The change of total spin 5 will qualitatively be the same 
as a change of temperature. Therefore we assume the 
spin to be constant: 5 = |. Furthermore we choose the 
zero point of energy to coincide with the band centers of 
the leads when no voltage is applied: BC me t = 0. Then 
there are seven remaining model parameters: the band 
widths of the metals W me t and the insulator Wi ns , the 
band center of the insulator BCi ns , the band occupation 
in the metals n, the exchange coupling J, the tunnel cou- 



pling e 



k A fkj 



and the number of insulator layers d. The 
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last two parameters play a special role because for the de- 
scription of a concrete experiment all the other ones are 
fixed. The number of layers can be determined approx- 
imately by the thickness of the sample and the lattice 
constant. The tunnel coupling is the only really free pa- 
rameter in the model and can be used to fit theoretical to 
experimental results. We assume that it is independent 
of the wave vectors eg ^ = e^f and choose the following 
ansatz for the tunnel coupling: 

e n L = e\ d e- x< - n -^ (nel,...,d) (20) 
e n R = e\ d e- x{d - n) (21) 

The exponential term is responsible for damping the cou- 
pling between the metals and insulator layers that are 
further away. X d will decrease the overall strength of the 
coupling with increasing thickness of the insulator for 
A < 1. The disadvantage of this ansatz is that one gets 
another undetermined parameter into the theory, namely 
A. One way to determine it is by comparison with other 
theories. Simmons^ derived a tunneling current formula 
through a thin insulating film. For intermediate voltages 
(0 < V < 0o, where 4>o is the average barrier height) it 
reads 

Ad) = |0-f )e- fl «- 

- (0o + ^) e - B^^ ^ / ^+*) (22) 

A and B are constants. By comparing the results of our 
theory with this formula, we find an almost perfect fit 
for A = 0.6. The dependence of the tunneling current on 
the remaining constant e is generally complicated, but it 
was found numerically that the current scales with e 4 for 
small values of e < 0.3 eV. Therefore it can formally be 
included in the unit of current together with the number 
of lattice sites in the insulator N^. 

[J] = ^ (23) 

In this section both leads consist of the same normal 
metal. Before beginning a systematic evaluation of the 
model it is useful to get a better understanding of how 
hybridization can be used to model tunneling. In Figs. 
[3] and d] the QDOS of the three regions are shown, with 
and without tunnel coupling e. Since there is no direct 
contact between the two leads, electrons can only get 
from the left to the right metal by hopping through the 
insulator. In the case e = this is not possible, because 
there are no allowed states in the energetic region of the 
occupied metal states. If the tunnel coupling is turned 
on, there will be allowed metal states in the energetic 
region of the insulator and vice versa. Therefore hopping 
between the layers becomes possible. Of course there still 
won't be a net current without an applied voltage. 
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FIG. 3: QDOS of a M/I/M system without tunnel coupling 
for T = 0. 
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FIG. 4: QDOS of a M/I/M-Systems with finite tunnel cou- 
pling for T = 0. The insulator states near the Fermi energy 
are additionally enlarged for clarity. 



A. One-layer insulator 

A typical current profile for a one-layer insulator is 
shown in Fig. El For small voltages the current increases 
in a linear way up to a maximum and later on it goes back 
to zero. Without an applied voltage there are only occu- 
pied or unoccupied states lying opposite to each other. 
Due to the Pauli principle there can't be any tunneling in 
this case. By increasing the voltage the right metal gets 
shifted downwards compared to the left one. For small 
voltages the number of occupied states on the left lying 
opposite of unoccupied states on the right grows more or 
less linearly. Therefore the current increases in a linear 
way, too. When there is a maximal number of occupied 
states lying opposite to unoccupied states the current 
reaches its maximum. For higher voltages it decreases 
again and will become zero when the voltage exceeds the 
band width of the metal (5 eV in this case), because then 
there is no overlap between the two bands anymore. This 
bevavior is the same for both spin directions, but obvi- 
ously the spin-up current is higher over the whole voltage 
range, i.e. there is a finite spin polarization. The reason 
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FIG. 5: Qualitative behavior of the spin-dependent tunnel- 
ing current for a one-layer insulator. Parameters: Wmet — 
5eV, Wins = leV, BC lns = SeV, n = 0.5, ^- = 0, J = 0.3eV 



behind this is the energetic distance between the metal 
and insulator bands. For the spin-up bands this distance 
is less than for spin-down, resulting in a stronger effect 
of the hybridization, which in turn will lead to a higher 
current. 

In Fig. [5] the dependence of the maximum of the to- 
tal current J = J* + on the model parameters is 
shown, again for a one-layer system. In the case of the 
exchange coupling and the temperature the insets also 
show the spin-dependent behavior of the maximum. For 
all the other parameters it is qualitatively the same as 
the shown curves. In the following list we will explain 
these curves: 

Band width of the metals W met : 
Increasing the width of the metal bands automatically 
decreases the maximum of the QDOS. This is shown in 
the inset. A decrease of QDOS means that there are 
less electrons which can tunnel and at the same time less 
empty states which they can tunnel into. Both effects 
reduce the current. 

Band width of the insulator Wins : 
By increasing the band width of the insulator and leav- 
ing its center of gravity constant the height of the tunnel 
barrier gets smaller. In the model this corresponds to a 
higher hybridization strength because the energetic dis- 
tance between metal and insulator bands decreases. So 
one expects an increase of current. 

Band center of the insulator BCi ns : 
The explanation of this effect is similar to the band 
width. The band center of the insulator is directly pro- 
portional to the height of the tunnel barrier. Shifting the 
band to higher energies will lower the current. 

Exchange coupling J: 
The dependence of the current on J can only be un- 
derstood by looking at the spin-dependent contributions 
shown in the inset. The spin-down current decreases 
slightly while the spin-up current shows a strong increase 
resulting in an increase of the total current. This behav- 
ior can be explained with the help of the SCj ns -curve. 



For J = eV both spin-down and spin- up bands have the 
same center of gravity. Increasing J shifts the spin-up 
band to lower and the spin-down band to higher energies 
(compare Figs. [7] and [SJ. It can immediately be seen 
from the picture that there will be a slight decrease in 
spin-down current and a bigger increase of spin-up cur- 
rent. 

Band occupation n: 
If the metal band is completely occupied (n=l for both 
spin directions) or completely empty, it becomes an in- 
sulator and the current vanishes. For small n there is 
only a small amount of electrons in the left metal which 
can tunnel. Vice versa, for high n there is only a small 
amount of unoccupied states in the right metal which 
can be tunneled into. For half filling there is an un- 
occupied state in the right metal for every electron in 
the left metal. Therefore one would expect a symmetric 
curve, with positive slope for n < 0.5, a maximum at 
n = 0.5 and negative slope for n > 0.5. But the result 
of the model calculation shows that the curve is shifted 
asymmetrically to higher band occupations. The reason 
behind this is that the distance between Fermi energy 
and the lower band edge of the insulator is smaller for 
high n. Therefore the symmetry between small and high 
n is broken. 

Temperature T/Tq: 
With increasing temperature spectral weight is shifted 
from the lower to the upper spin-up band (compare Fig. 
in]) . This effectivly increases the tunnel barrier for spin- up 
electrons and thus the current is decreased. With a sim- 
ilar argument it can be understood why the spin-down 
current will increase with increasing temperature. At 
the Curie temperature Tc both curves meet and above 
Tc the current is practically no longer dependent on tem- 
perature. 

The explanation of the J-curve with Fig. [7| will 
be crucial for the understanding of the spin polarization 
P = jrzji ■ Since the slope of the BCi ns -cuvve decreases 
drastically for higher band centers, it follows that higher 
band centers (or tunnel barriers) are bad from the point 
of view of the spin polarization. 

In Fig. [TO] the spin polarizations are shown as functions 
of the applied voltage and the model parameters. The 
spin polarization decreases with increasing voltage for all 
parameters. In the case of one insulating layer only the 
band center of the right metal is influenced by the volt- 
age, while the band centers of the insulating layer and 
the left metal remain constant. The voltage shifts the 
right band center to lower energies, i.e. the height of 
the tunnel barrier on the right side increases. With the 
argument given in the paragraph above this leads to a 
lower spin polarization. The dependence on the model 
parameters will be discussed in the following list: 

Band width of the metals W me t- 
The spin polarization only shows a weak dependence on 
W met - It changes only about 2 % during a change of 
band width from 3 to 15 eV. The apparent sudden tilt 
for high voltages and low band widths results from nu- 
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FIG. 6: Dependence of the tunneling current maximum on the model parameters. The inset of the VFmet-figure shows the 
maximum of the metallic QDOS in dependence on the band width. The insets of the 7^7- and J-figures show the contributions 
from the two spin directions. For all the other parameters they behave like the total current and are not shown for clarity. 
The following standard parameters were used: Wmet = 5eV, Wi ns ~ leV, BCi na = 3eV, n = 0.5, T = OK, J — 0.3eV. For each 
picture just one of them was varied. 
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FIG. 7: (Color Online) Graphical explanation of the behav- 
ior of the current maximum in dependence on the exchange 
coupling J with the BCins-curve from Fig. [6] 
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FIG. 8: (Color Online) Insulator QDOS for different values 
of the exchange coupling constant J without coupling to the 
metals. Parameters: Wi ns = leV, BCi„ s = 3eV, S = |,T = 
OK 



merical difficulties since the current is almost zero for 
those parameters. 

Band width of the insulator Wi ns : 
High insulator band widths are good for spin polariza- 
tion. This is especially obvious for small voltages. For 



higher voltages this increase is not very pronounced. 

Band center of the insulator BCi ns : 
The spin polarization changes drastically in dependence 
on BCi ns i especially for the case of low band widths (no- 
tice the spin polarization scale!). For low voltages it 
grows by a factor of three from about 20 % to almost 
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FIG. 9: (Color Online) Insulator QDOS for different values of 
temperature t^j without coupling to the metals. Parameters: 

Wins = leV, BC ins = 3eV, S = | , J = O.SeV 

60% when decreasing the band center from 5 to 1 eV. 

Temperature T/Tc'- 
The spin polarization closely resembles the behavior of 
the magnetization, i.e. it has a saturation value for 
T = 0, slowly decreases with increasing T and falls off 
very rapidly near the Curie temperature. Above Tc it is 
zero, as expected. 

Band occupation n: 
The spin polarization is constant for low n. For very high 
n and low voltages there is a clear increase since the tun- 
nel barrier is lowered with increasing n. This is good for 
high polarization as explained above. 

Exchange coupling J: 
The spin polarization shows the strongest change in de- 
pendence on J. By increasing the exchange coupling one 
can reach polarizations of almost 80 % (although such a 
high J might not be realistic with the chosen band width 
of 1 eV). The increase is approximately linear in J. 



B. Two-layer insulator 

If we replace the one-layer insulator by a two-layer in- 
sulator the dependence of the spin polarization on the 
model parameters remains the same but the dependence 
on the voltage is reversed. This is shown exemplary in 
Fig. [II] for BCi ns . Now the spin polarization increases 
with increasing voltage. If we denote the metals by L 
and R as above and the left insulator layer by II and the 
right layer by 12 the following happens when the voltage 
is increased: the distance between L and II and R and 
12 stays the same. At the same time the distance be- 
tween L and 12 is decreased and the distance between R 
and II is increased. The first effect will increase the spin 
polarization (lower barrier height!) while the second one 
will decrease it. Due to the high slope of the BCi ns -curve 
for low band centers the first effect dominates and leads 
to the observed behavior of the spin polarization. 



Insulators with more than two layers behave qualitatively 
in the same way as discussed here. Only the strength of 
the current falls off with increasing number of layers. 
DeWeert and Girvin 20 performed a Boltzmann-equation 
study of M/I/M-tunneling where they also found that 
spin polarization will increase with voltage. This is 
not surprising since they made assumptions similar to 
ours. However, they treated the ferromagnetism in a 
phenomenological way by introducing different tunnel- 
ing life times for the two spin directions. In our case 
the magnetism is a direct consequence of the considered 
microscopic model. 



C. Comparison to experiment 

Now we want to compare our theoretical results with 
experimental measurements for an Al/EuS/Al-system 2 ^. 
In Fig. [T2]the results of the experiment (circles and tri- 
angles) are shown together with our calculation (lines). 
The model parameters are taken from Wachter— . The 
lattice constant of bulk EuS is ao = 5.97 A. In the exper- 
iment the insulator had a thickness between 17.6 A and 
21.9 A which amounts to approximately three layers. Of 
course this is an approximation since the lattice constant 
will change for small thicknesses, but the order of mag- 
nitude should be the same. The parameters for Al were 
taken from Papaconstantopoulos 23 . The tunnel coupling 
e is the last unknown parameter. The unit of current is 
jre i Ni, therefore e can be used to fit the intensity of the 
theoretical current to the experimental results. This was 
already done in Fig. QJwith e 4 N, = 5.3- 10 6 . The surface 
area of the sample was A = 3 • 10~ 3 cm 2 . Assuming that 
the lattice sites are ordered quadratically one can esti- 
mate the number of lattice sites to Ni = 8.4- 10 11 . Hence 
e is known: e w 0.05 eV. Of course, this is a strongly 
simplified approximation, but it shows that the tunnel 
coupling is small enough to justify its inclusion in the 
unit of current, which is only possible for e < 0.3 eV. 
There are several possible reasons for the quantitative 
difference between theory and experiment. Aluminum is 
not well described by a single band tight-binding approx- 
imation. Since the metal QDOS has a strong influence on 
the current characteristics one has to expect differences. 
Another reason is the treatment of the tunnel barrier. 
It was assumed to be rectangular and was rigidly tilted 
by the applied voltage. This assumption is, especially 
for reasons of simplicity, often made. But on the other 
hand it is known to be a very crude approximation for 
real barriers since they can have varying barrier heights 
and even holes 2 ^. Another quantitative difference is the 
behavior in dependence on temperature. Both in exper- 
iment and theory the slope of the curves decreases with 
increasing temperature, but this decrease is stronger for 
the experimental curves. The origin of this difference 
most likely lies in the treatment of the tunnel coupling. 
It was already shown that spin-up and spin-down current 
approach each other with increasing temperature, but 
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FIG. 10: (Color Online) Spin polarization as a function of voltage and model parameters. Blue (dark gray) corresponds to 
low, yellow (white) to intermediate and red (light gray) to high polarization. For better clarity it is marked by arrows in which 
direction the scales increase. Parameters: W m et = 5eV, Wins = leV, BCi ns = 3eV, n = 0.5, T = OK, J = 0.3eV 



that the spin-up current decreases a little bit faster than dence of the tunnel coupling which was neglected in this 

the spin-down current increases. Thus the total current paper. 

decreases, but not enough to explain the experimental 

results. This decrease would be larger if the BCi ns curve 

in Fig. would be steeper. That might be accomplished 

by explicitly taking into account the wave- vector depen- 
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FIG. 11: (Color Online) Spin polarization for a two-layered 
insulator as a function of voltage and the band center of the 
insulator BCins- The arrows mark in which direction the 
scales increase. Parameters: W me t = 5eV, Wins = leV, n — 
0.5, T = OK. J = 0.3eV 
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FIG. 12: Comparison of the theoretical results (lines) with 
experimental measurements by Hao et al.— (circles and tri- 
angles). Parameter for the calculation: Wmet = 26eV, Wins = 
0.9eV, BC ina = -1.3eV, n = 0.33, J = O.leV, d = 3 



IV. METAL/INSULATOR/SUPERCONDUCTOR- 
SYSTEM 



In this section we replace one of the metals with a 
superconductor (M/I/M — * M/I/S). For the derivation of 
the tunneling current formula we assumed that the leads 
are non-interacting. Obviously, this assumption does not 
hold in the case treated here. Hl has to be replaced by 
the BCS-Hamiltonian 2 ^: 



H T . -» Hr 



C k M er Ck iW cr 



MeL.R 



A 2 

V 



A is the bandgap of the superconductor. The problem 
with this approach lies in the fact that Hl does not com- 
mute with the particle number operator anymore. In 
principle one has to derive a new current formula with 
this new Hl as was done e.g. by Zeng et al4. Here we 
use an alternative approach, which is known as the ef- 
fective medium approach. The basic idea is to split off 
the interacting part and use it for a renormalization of 
the one-particle energies and therefore for the density of 
states: 



(S)l 



(24) 



p( s \E) is the superconducting QDOS. Thus the model 
Hamiltonian remains the same, only the QDOS has to be 
changed in the numerical treatment. The disadvantage of 
this approach is the restriction to single electron effects, 
i.e. the neglicence of Cooper pair tunneling for example, 
but the comparison with experiment will justify this. 
In recent years a lot of research was focused on the trans- 
port properties of ferromagnet-superconductor junctions. 
Especially the role of Andreev reflection has been inves- 
tigated by many author s 2627282930 . Since Andreev reflec- 
tion involves more than one electron it cannot be mod- 
eled by our effective medium approach. Again this can 
be justified by comparison with experiment later in this 
section. 

The main difference between normal metals and super- 
conductors is the appearance of the bandgap. Its width 
is usually of the order of meV. Therefore we can restrict 
the discussion to voltages which correspond to these en- 
ergies. For higher voltages the results will coincide with 
the M/I/M curves. The reason for introducing the multi- 
layer insulator was the modelling of the band tilting by 
the applied voltage. In the case considered here this 
should have no effect on the current, since the tilt is of 
the order of meV with a total barrier height in the or- 
der of eV. The change from a one-layer to a many-layer 
system will merely change the intensity but not the qual- 
itative behavior of the current. 

For better comparison with experiment later on we fo- 
cussed on the differential conductance Hrr rather than 

dV 

the current itself. It is normalized on the total conduc- 
tance of the M/I/M system 



dj" 



dVMIS 

dJT+dJl 

dVMIM 



(25) 



In Fig. [13] the conductance is shown for several values of 
the Debye temperature Tp. For the parameters used we 
get the following bandgaps: 



A(T D = 15K) = 0, 46 meV 
A(T D = 30K) = 0, 92 meV 
A (To = 60K) = 1,85 meV 



(26) 



There is no current for voltages lower than these 
bandgaps. Furthermore the conductance curves are not 
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FIG. 13: (Color Online) Differential conductance of the 
M/I/S system for different values of the Debye temperature. 
Dotted lines are spin-down, continuous lines are spin-up con- 
ductances. All curves are normalized according to Eq. (125 p . 
Parameters: W me t = 3eV, W^ ns = leV, BCj ns = 3eV, e — 
0.03eV, T = OK, n = 0.5, J = Q.ieV 



symmetric about the origin: for positive voltages the 
spin-up current starts first followed by the spin-down cur- 
rent, for negative voltages vice versa. Ignoring this spin 
splitting the basic form of the conductance can be under- 
stood with the following arguments: for V = the chemi- 
cal potentials on both sides are equal and current cannot 
flow. With increasing voltage there still won't be any 
current since there are no states available in the super- 
conductor. Only when the voltage exceeds the bandgap 
\eV\ > A there will be occupied and unoccupied states 
opposite to each other. Hence one expects an increase 
of conductance at the edge of the bandgap. Since there 
are many available states above it this increase will be 
stronger than in the M/I/M system. 

The spin splitting of the conductance curves is due to 
the coupling between the ferromagnetic insulator and the 
non-magnetic superconductor. The hybridization repels 
the two bands. The spin- up insulator band is closer to 
the superconductor band and therefore causes a larger 
band shift. This can be seen from the QDOS in Fig. [TJ] 
For positive energies the band edge of the spin-up band 
is just below the bandgap while the spin-down band is 
just above it. When applying a positive voltage one ex- 
pects the spin-up current to start first which was indeed 
observed in Fig. [T3] The QDOS of the insulator is finite 
inside the bandgap, i.e. there will be a finite contribution 
to the tunneling current from this region! This effect has 
its origin in the hybridization with the (constant) QDOS 
of the metal and is physically not explainable. Though 
the current inside the bandgap is finite, it is still at least 
three orders of magnitude lower than the current above. 
Thus it cannot be seen in Fig. [T51 

Compared to the M/I/M system the spin polarization 
shows a richer structure (Fig. fTS"]) . Due to the different 




i i i I i i Hi i i 
2.5 2 1.5 1 0.5 

QDOS (10 3 ) QDOS (10- 7 ) 

FIG. 14: (Color Online) QDOS of the superconductor and 
the insulator in the range of the bandgap. There is a clear 
splitting between spin-up (green/light gray) and spin-down 
(black) observable. Parameters as in Fig. [13] 
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FIG. 15: (Color Online) Spin polarization of the M/I/S 
junction for different values of the Debye temperature. The 
dotted region around the origin was omitted since the spin 
polarization is not defined there. The same parameters as in 
Fig. [13] were used. 



onsets of the current for the two spin directions there 
are two strong peaks at ±A. For positive voltages one 
observes positive spin polarization with a maximum of 
about 90 %. It decreases slightly with increasing width 
of the bandgap. For negative voltages the spin polariza- 
tion will be negative in a narrow region with a minimum 
of about -50 %. 

The tunnel coupling e determines the strength of the 
hybridization between the leads and the insulator. The 
bigger it gets, the stronger the spin splitting of the su- 
perconductor should be. This effect is confirmed in Fig. 
[TBI For e — 10 meV there is almost no spin splitting 
observable and the spin polarization only has two small 
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FIG. 16: Differential conductance of the M/I/S system for dif- 
ferent values of the tunnel coupling e. Dotted lines represent 
spin-down, continuous lines represent spin-up conductances. 
They are normalized according to Eq. (|25|l . Parameters: 
WVet = 3eV;W ins = W,BC ins = 3eV,T D = 30 A", T = 
0K,n = 0.5, J = 0.3eV 
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FIG. 17: (Color Online) Spin polarization of the M/I/S 
system for different values of the tunnel coupling e. The spin 
polarization is not well defined in the dotted region around 
the origin. The same parameters as in Fig. [TS]were used. 



bumps accordingly. By increasing e the spin splitting is 
increased until the bandgaps of spin-up and spin-down 
almost do not coincide anymore for e = 100 meV. Thus 
in this case the highest spin polarization of about 95 % is 
reached which is also stable across a broad voltage range 
of about 1 meV. 

The last model parameter which will be discussed is 
the temperature T. In Figs. [TBI and [T9l two conductance 
curves and the according spin polarizations are shown for 
T = and T > 0. Increasing the temperature has two 
effects: first the bandgap is reduced which is equivalent 




-i o i 
Voltage (mV) 

FIG. 18: (Color Online) Differential conductance of the 
M/I/S system as a function of the applied voltage for two 
different temperatures. Continuous lines represent spin-up, 
broken lines spin-down conductances. Parameters: W me t = 



3eV, W h 
0.3eV 



leV,BC ; - 



3eV,Tn = 30K,n = 0.5, J 
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FIG. 19: (Color Online) Spin polarization of the M/I/S junc- 
tion for two different temperatures. The spin polarization is 
not well defined in the dotted region around the origin. Pa- 
rameters as in the figure above. 



to a reduction of To and thus nothing new. Second the 
Fermi edge is softened and this can lead to finite occupa- 
tion above the bandgap. Therefore current can already 
flow for voltages below the bandgap which leads to the 
observed softening of the conductance curves. For the 
same reason the maximum of the conductance is lower. 

The results of the model calculation will now be com- 
pared to experimental results obtained in the group of 
Mooder a 213132 . They used several europium chalko- 
genides as tunnel barriers. All the results are qualitativly 
the same, thus we will compare our results with the EuS 
measurements^ only. The tunnel junction consisted of 
an Au/EuS/Al system where Al becomes superconduct- 
ing below T — 1.2 K. In Fig. [^D] the experimentally 
measured conductance is shown for several values of the 
applied magnetic field. Even without applied field they 
measured a finite Zeeman splitting which corresponds to 
an internal field of B* = 0.5 T. The dependence of this 
internal field on the applied field is shown in the inset 
of the figure. Obviously the internal field is much bigger 
than can be explained by the applied field alone. This is 
a strong indication of the coupling between EuS and the 
Al-lead. 
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FIG. 20: Experimental conductance curves for an Al/EuS/Au 
tunnel junction and different strengths of the applied mag- 
netic field (from Hao et al. 21 ). The inset shows the depen- 
dence of the internal magnetic field B* on the applied field 
H. 
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FIG. 22: Contributions to the conductance from two differ- 
ent domains of the EuS (from Hao et alM-). (a) shows the 
schematic design of the Al/EuS/Au system, (b) the Zeeman- 
split density of states of the Al and (c) the expected conduc- 
tance. Its overall behavior is the same for both domains, only 
the spin directions are exchanged. 





t 


1 


1 1 








♦ 1 




1 1 








i V — 


, J 


1,1, 



-3 -2-10 1 

Voltage (mV) 



density of states according to Fig. [22] Obviously there is 
a good qualitative agreement between their expectation 
and the results of our model calculation. 
Andreev reflections would show in the conductance curve 
as a characteristic peak near V = (compare e.g. the 
conductance curves measured by Soulen et ali^I). This 
peak was not measured by Hao et al. and it cannot ap- 
pear in our model as was explained above. Thus we can 
draw the conclusion that Andreev reflection is strongly 
suppressed in this system which is a consequence of the 
high spin polarization. Hence only single-particle exci- 
tations contribute to the conductance and the effective 
medium approach should be sufficient. 



SUMMARY AND OUTLOOK 



FIG. 21: Normalized total conductance as a function of the 
applied voltage. The arrows mark the spin direction which 
contributes most to the respective maximum. The width 
of the bandgap is about 0.9 meV. Parameters: W me t = 
3eV,Wins = leV,BC ins = 3eV,n = 0.5, T = 0, e = 
70meV, J = 0.3eV, T D = 30^ 



In our model the internal field is created by the hy- 
bridization between the ferromagnetic insulator and the 
paramagnetic leads. Therefore an increase of the tunnel 
coupling in the model corresponds to applying an exter- 
nal field in the experiment. In Fig. [21] the theoretical 
total conductance is shown for e = 70 meV. The model 
calculation for the smaller e = 30 meV in Fig. [TH] agrees 
qualitatively with the H = T-curve. Furthermore there 
is also an agreement between the curves with applied ex- 
ternal field and the results for e = 70 meV. The authors 
also discuss their results in terms of the superconductor 



The topic of this paper was the theoretical modelling of 
spin filter experiments, i.e. the spin-dependent tunneling 
through a ferromagnetic insulator. Due to the coupling 
between two outer leads and the insulator with in general 
different chemical potentials the tunnel junction is not in 
thermal equilibrium anymore. Thus one has to use the 
non-equilibrium extension of the usual many-body the- 
ory, the so-called Keldysh formalism. We presented a 
model which can be used to simulate tunneling of elec- 
trons through a potential barrier in the presence of in- 
teractions. The insulating region was described by the 
Kondo Lattice Model and treated within an interpolat- 
ing self-energy approach. It can be modelled as a one- 
or multi-layered system. The leads consisted of either 
two normal metals or one normal metal and a supercon- 
ductor. The metals were assumed to be non-interacting, 
while we used the BCS theory to describe the supercon- 
ductor. The tunneling process itself was simulated by a 
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FIG. 23: Time contour C in the Keldysh formalism 

hybridization between the conduction bands of the leads 
and the insulator. Following Haug and Jauho^ we de- 
rived a current formula which holds for non-interacting 
leads and arbitrary interactions in the insulator. In all 
considered systems a finite spin polarization was found. 
The degree of polarization was strongly dependent on the 
model parameters. It was possible to explain the current 
characteristics with means of the QDOS. The comparison 
with experiment has shown that the model predicts the 
right qualitative behavior, but that there are also clear 
quantitative disagreements. They stem from the not suf- 
ficiently realistic treatment of the tunnel coupling and 
the band structure of the used materials. It was espe- 
cially suspected that the tunneling probability decreases 
too slowly with increasing barrier height. Therefore the 
spin polarization is probably underestimated. 
The superconducting contacts were treated as an effective 
medium. A strong spin splitting of the superconductor 
QDOS was observed which was due to the coupling with 
the ferromagnetic insulator. Thereby obvious maxima 
and minima of spin polarization formed near the edges 
of the bandgap. Even negative polarization could be ob- 
served. Comparison with experiment showed a qualita- 
tive agreement with the model results for the M/I/S sys- 
tem, too. 

There are several aspects of the model which might be 
improved. One is the treatment of the tunnel coupling. 
Especially by the comparison with experiment for the 
M/I/M system it became clear that the assumption of 
wave-vector independent tunnel coupling is too much 
simplified. It would be very desirable to have a method 
to derive the tunnel coupling from the other model pa- 
rameters which are easier accessible by experiments. A 
second important improvement would be the use of real 
bandstructures, e.g. from DFT calculations. 

I 



Putting this expression in (|A1[) and Fourier transforming 
the whole equation leads to the result Eq. ((9]). 



APPENDIX A: DERIVATION OF THE 
CURRENT FORMULA 

Evaluation of the commutator in Eq. ([5]) with H given 
by Eqs. ©-(SI) leads to 

JiW = f E Re (Ck^« ff W^a,n(«))) 

= | E K^ kj G< w (M)) (Al) 

Following Keldysh 4 wc defined the lesser Green's function 
in the last step: 

Gw (Tl (i,0 = i(ci L(T (t')d kl(Tin (t)) (A2) 

Obviously the current is known if it is possible to cal- 
culate this Green's function. This task can be ac- 
complished by deriving the equation of motion of the 
corresponding causal Green's function G kjkiCT n (t, t') = 
—i{Tdk I tj,n(i)c^ cr (t')). Its formal solution is 

oo 

Gkjkitr.n (*)*') = t / dt l E £ k*ki " 

— oo 

• G^ J n m (*.*i)flk i (*i>0 (A3) 
where we have used the free lead Green's function 

(-ift^7 - ^) & L (*, = H6(t - t>). (A4) 

Since the non-equilibrium Green's function (NEGF) has 
formally the same perturbation expansion, with integra- 
tion over the real time axis changed to integration over a 
time contour C shown in Fig. 1231 ^, we can immediately 
write it down as 

G f k J k£<r,n(*)* / ) = ^ J dr E e k*kt ■ 

c 

■ Gk I a,nm{t,T)gk L {T,t') (A5) 

Applying the analytic continuation rules given by 
Langreth 1 ^ to get the lesser part of the NEGF yields 



(A6) 

I 

APPENDIX B: CALCULATION OF THE 
INSULATOR GREEN FUNCTIONS 



OO 

= \ I *i E (GUnm(*. h)g< L (t u t>) + G$ mnm (t, h)gZ L (t x , f )) 



The calculation of the retarded and lesser Green's func- 
tion in the insulator runs along similar lines as the deriva- 
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tion of the current formula. The first step will be to 
solve the equation of motion of the corresponding causal 
Green's function. The Hamiltonian is again given by Eqs. 



dU)- <j4j) . This time we have to explicitly take the inter- 
actions, i.e. the KLM, into account. Its self-energy was 
defined by Eq. (fT3j) . Then the equation of motion reads 



Mk/^kA/k/cr.m 

(E) (Bl) 

l l M k M 

M£L,R 



The 'mixed' Green's function G£ Afkj . CT m (E) can be cal- The equation of motion becomes 
culated as in (|A3|) . The final expression is very similar: 

G k M kia,m(E) = ■^X! e kMk ^ 5 , k L ( ^; ) G 'k /^ T,^ m ( £, ) ( B2 ) 

i 



E (ESm e£ - SZ^raW) G^ lm {E) = H5 nm + ^ E £ ^UAm^ {E)G^ lm (E) (B3) 



l " M k M / 



On the l.h.s. we defined the free insulator Green's func- 
tion. One should note that it includes the interactions 
through the ISA self-energy. So 'free' does not mean 
that it does not contain interactions but that it solves the 
equation of motion for vanishing tunnel coupling. There- 
fore it can be calculated in thermal equilibrium. This is 
the essential premise for using the ISA self-energy within 
this formalism. 

The equation of motion can be brought into the form 
of a Dyson equation by defining a tunnel self-energy as 
follows 



The next steps of the calculation are to Fourier trans- 
form this equation to time dependence, replace integra- 
tion over the real time axis by integration over the time 
contour and then to use the analytic continuation rules 
again to get the lesser and retarded part of the NEGF. 
To keep the notation as simple as possible we will use 
a matrix notation from now on where matrices will be 
marked by a hat. Then the Dyson equation for the re- 
tarded Green's function reads: 



A,, 



i{E) 



1 



£ £■ 

M k M 

MeL,R 



k»/k 



^9i M {E) (B4) 



With this definition the equation of motion reads 

Guja^nmiE) = <?kj a,nm (E) + 

+ E 9Zi<r,np( E )^ii ,pq(E)G^ ia qm (E) 

pq 



Gl ia {E) 



9lAE) +g r kl AE)K 1 (E)Gl irT (E) (B5) 
Similarly the lesser Green's function solves the Keldysh 
equation: 



G£ ia (E) = (i + G kia (E)A ki (E))g< a (E) (i + K i (E)G a kl (E))+Gi ia (E)A< i (E)G a kia (E) (B6) 

I 

The free lesser insulator Green's function g kir7 (E) can non-equilibrium, but as was pointed out before, the free 
be expressed by the free retarded Green's function with Green's functions can be calculated in equilibrium. One 
the help of the spectral theorem. It doesn't hold in 
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& a (E) = -2if ins (E)Img r kia (E) (B7) 
where fins(E) is the Fermi function in the insulator. 
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